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The variance in the winding number of various random fractal curves, including the self-avoiding 
walk, the loop-erased random walk, contours of FK clusters, and stochastic Loewner evolution, have 
been studied by numerous researchers. Usually the focus has been on the winding at the endpoints. 
We measure the variance in winding number at typical points along the curve. More generally, 
we study the winding at points where k strands come together, and some adjacent strands may be 
conditioned not to hit each other. The measured values are consistent with an interesting conjecture. 

PACS numbers: 05.50.+q, 64.60.Fr, 64.60.Ak 



Duplantier and Saleur [l[ studied the winding angle 
between the two endpoints of a finite self-avoiding walk 
(SAW) in 2D, and indeed, a broader class of curves. Us- 
ing exact but nonrigorous Coulomb gas methods, they 
found that the distribution of winding angle approaches 
a Gaussian and they explicitly computed the variance. 
When the endpoints of the walk are distance L apart, 
the winding variance is ~ (8/g) logL, where g is a model- 
dependent parameter which is 3/2 for SAW. The winding 
angle at a single endpoint (relative to the global average 
direction of the curve — see below for a precise definition) 
is a Gaussian with variance (4/g) logL We found ex- 
perimentally that the variance in the winding at typical 
(random) points along the curve was only 1/4 as large 
as the variance in the winding at the endpoints. More 
generally, when k strands of the curve come together at 
a point, the winding angle variance is 1/fc 2 as large as at 
the endpoints; Eq. Q below generalizes this further. 

Remark: After our initial experiments we learned that 
the 4/(<?fc 2 )logL formula is also contained in unpub- 
lished notes of Duplantier. However, to our knowledge 
the winding at typical points or points where k strands 
come together is not mentioned in the literature, except 
in the case of loop-erased random walk (LERW), where 
we have identified a minor oversight in the calculations 
that resulted in incorrect values being reported. Our ex- 
periments can be seen as a test of Duplantier's Coulomb 
gas predictions. We also report on other random fractal 
curves, for which the Coulomb gas methods do not apply. 

The main object of our study is the 2D Fortuin- 
Kasteleyn (FK) random cluster model @, H| at critical- 
ity, specifically the contours of the clusters. The FK 
model is like bond percolation with edge probability p, 
but there is another parameter q, and Pr[configuration] oc 

p # bonds ^ missing bonds ? # clusters^ For each ^ there 

is a critical p above which the system percolates. When 
P = ^critical , the contours of the clusters form a system of 
loops called the fully- packed loop model (FPL) [4H6(. See 
Figure 1, left panels. A loop configuration occurs with 
probability proportional to loops where n — Jq ■ 
In addition to the perimeter (or "hull") of a cluster, 
the "external perimeter" has also been studied (see Fig- 
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FIG. 1: The contours and external perimeter of critical FK 
clusters. The upper left panel shows a FK random cluster 
configuration with q — 3 at criticality. On the upper right 
panel we closed off the narrow passageways of this FK config- 
uration by connecting adjacent pairs of vertices that belong to 
the same connected component. In the lower panels we show 
the fully packed loop configurations that come from the above 
bond configurations. The loops on the left traverse the hull 
or perimeter of the clusters. The loops on the right traverse 
the external perimeter of the clusters. The longest perimeter 
loop and longest external perimeter loop are shown in bold. 



ure 1, right panels). Grossman and Aharony 7] found 
experimentally that by closing off narrow passageways 
on the hull of percolation clusters, the fractal dimension 
drops from 7/4 0] to 4/3, and that furthermore the pre- 
cise definition of "narrow" had little or no effect on the 
fractal dimension 4/3. See Q for an explanation of this 
phenomenon in terms of path crossing exponents. The 
external perimeter of percolation clusters is believed to 
be essentially the self-avoiding walk 0,0- U|, and the ex- 
ternal perimeters of FK clusters for other values of q are 
also interesting. Therefore in addition to studying the 
hulls of the FK clusters, we also studied their external 
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q: cluster fugacity in FK random cluster model n 2 = q [4-^6j (perimeter only, not external perimeter) 

n: loop fugacity in fully packed loop model n = — 2cos(7rgi) [12. 131 

g: coupling constant of associated Coulomb gas « = K\ = 4/gr [lj 

Kfc: winding angle variance when k curves meet at a point g 

ext. perimeter 

= 1/9 

perimeter 

[14] 

Df. fractal dimension of curves Df = 1 + l/(2g) = 1 + k/8 [12. 151 



TABLE I: Summary of Kk measurements for the perimeter and external perimeter of FK clusters. The table goes up to 
q — 4, beyond which the FK model has a first order phase transition without large loops (l6l . Il7| |. The values of the fractal 
dimension D / are summarized from [?rM Q H-H^] , and Ki comes from [l[ . The values for K2 for the perimeter and external 
perimeter were measured as described below. K3 for LERW was measured by looking at the "triple points" of uniformly random 
spanning trees - points where three strands come together. Ki and «6 for the spanning tree perimeter are K2 and K3 for LERW 
respectively. For q = 1, K3 is ki of the external perimeter @, and likewise K5 = n'4 (defined below). The measurements are 
consistent with the hypothesis that Kk — ni/k 2 , a formula that also appears in unpublished notes of Duplantier [25| . 



perimeters by closing off passageways of lattice spacing 
1 and looking at the hulls of the resulting clusters. 

The perimeters and external perimeters of FK clusters 
are closely related to a variety of models, most notably 
the stochastic Loewner evolution (SLE) process intro- 
duced by Schramm [2^]. In a discretized version of SLE, 
a curve in the plane grows as follows: the portion of the 
plane not in the curve is conformally mapped to the half 
plane, with the tip of the curve mapped to the origin. 
A small random cut is then made starting at the origin, 
where the slope of the cut is controlled by a parameter 
k. The original curve gets extended by the pre-image of 
this small cut, and the process repeats. The variance in 
the winding angle at the endpoint of SLE K is k log L [3] . 

The SLE process describes the limiting behavior of a 
variety of statistical mechanical models in 2D. Schramm 
proved that SLE2 gives the scaling limit of loop-erased 
random walk (LERW), provided that LERW has a 
conformally invariant scaling limit; recently Lawler, 
Schramm, and Werner [53] proved this without assump- 
tions. Smirnov [28| proved that critical site percola- 
tion in the triangular lattice converges to SLEg. Lawler, 
Schramm, and Werner [29[ proved that the frontier of 
Brownian motion (with suitable boundary conditions) 
converges to SLEs/3. There are good theoretical [30( 
and experimental [31] reasons to believe that SLE 8 / 3 also 
describes the self-avoiding walk (SAW) . Calculations by 
Kenyon and Schramm [32j suggest that SLE4 describes 
the loops arising from superimposing two domino tilings. 
Rohde and Schramm [l5j proved that the fractal dimen- 
sion Df of SLE K is at most l + re/8 when k < 8, and their 
calculations suggest Df = 1 + k/8. See also 18, 2(J 33- 



37[ for further results on SLE. Schramm conjectured that 
the contours of FK clusters at criticality have the same 
local properties as SLE K , where k depends on q. 

We study the perimeter and external perimeter of FK 
clusters by looking at the winding angle function. Given 
a loop or a path in the plane or on a torus, we define the 
winding angle function w(), a function of the edges, as 
follows. We pick an arbitrary starting edge e on the loop 
or path and an arbitrary value for the winding function 
w(e) at that edge. The winding at a neighboring edge e' 
is defined by w(e') = w(e)+ the turning angle from e to 
e' measured in radians. This definition applies to paths 
or noncontractable loops, i.e. loops that wind around the 
torus. If a loop is contractable to a point, this would 
yield a multivalued winding function, so we adjust the 
definition of w(e') — w(e) by 27r/ (length of loop) to get a 
single-valued winding function. This specifies the wind- 
ing angle function up to a global additive constant; we 
choose the value of this global constant to make the av- 
erage winding angle of the edges on the loop or path 0. 

When k strands of the perimeter or external perimeter 
converge on a point, the winding angle variance should 
scale clS Kfc log L. Table U summarizes our measurements 
of Kk, suggesting Kk = Ki/k 2 (see also 2^, 38]). 

Remarks on LERW: Our simulation values for K2 and 
K3 for LERW disagree with the values previously re- 
ported by Kenyon [39( by a factor of 4 and 9 respectively. 
These calculations used the Temperley [4(J correspon- 
dence between spanning trees and dimer systems, and 
Kenyon correctly and rigorously computed the variance 
in the height function of the associated dimer system 
when there were 1, 2, or 3 paths approaching a point. 
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The height function of the dimer system is related to the 
winding angle for the paths: when there are k paths, 
each winding changes the dimer height function by 4fc. 
The factor of k was omitted, leading to the factor of k 2 
discrepancy in the winding angle variance. 

The longest contour of an FK configuration is likely 
to hit itself many times (which is why the perimeter 
and external perimeter arc different); the places where 
the contour hits itself are called pinch points. For the 
longest contour we identified the pinch point giving rise 
to the longest pinch. At this point there are four strands 
that travel a distance on the order of the box length L, 
suggesting that the winding angle variance at this point 
should grow as K^logL. However, as noted by Schramm 
24] , the pinch point with the longest pinch is an atypical 
pinch point because there are two adjacent strands are 
conditioned not to hit each other — if they did hit each 
other, then this would create a longer pinch. Thus the 
winding angle variance at the longest pinch point is gov- 
erned by a different constant k 4 , and grows as K 4 log\L. 

In general let n' k be the winding angle variance coeffi- 
cient when there are k strands meeting at a point and two 
adjacent strands do not hit each other (when k = 2, the 
left side of one strand may hit the right side of the other 
strand, but not vice versa). When q = 4, the strands do 
not hit each other anyway [15| . so n' k = Kk- When q = 1 
and k strands meet at a point, conditioning two adjacent 
strands not to touch has the same effect as adding an ex- 
tra strand: K k — Kk+i @, EH- For other values of q it is 
plausible that requiring two of the strands not to hit each 
other has the effect of adding some fractional number of 
strands f(q) between the strands required not to hit each 
other; similar phenomena have been observed elsewhere. 
When two strands meet at a point that happens to be 
on the external perimeter, the right side of one strand 
does not hit the left side of the other strand. Thus k' 2 
for the perimeter is K2_for the external perimeter, which 

(by K 2 = Kl/4 = I/3 [H and g ext . perimeter = 1/Sperimeter 

[14j ) in turn is 1/^2 (for the perimeter), giving 

Kl /(2 + f(q)) 2 =4//Si 

f{q) = K 1 /2-2 

4 =Kl/(fc + K 1 /2-2) 2 . 

For example, when q = this predicts k' 4 — 2/9. Indeed, 
the largest pinch point for SLE§ corresponds to a triple 
point of the spanning tree, for which we already have the 
value 2/9. For other values of g, our measured values of 
K4 appear to be consistent with this formula. 

More generally, when k strands meet at a point, and 
j adjacent pairs do not hit each other, we expect the 
winding angle variance to grow like 



logL. 



(*) 



multiple of 4, 5, 6, or 7, starting with L = 4x2 3 ,5x2 3 ,6x 
2 3 ,7 x 2 3 ,4 x 2 4 ,...), we generated 10000 random FK 
configurations using the methods described in [4^, |42| |. 
In each one we identified the longest contour (and also 
the longest outer contour) and computed the winding 
angle function as defined above. For «2, the square of 
the winding at a random single edge on the loop is an 
estimator of the winding angle variance of the loop, but 
a more efficient estimator is the average square of the 
winding angle of edges on the loop. For « 4 , we measured 
the square of the winding at the pinch point of the longest 
pinch of the longest contour. The data for the longest 
contour when q — 1 (percolation) is shown in Figure [5J 
whose caption explains how we estimated /t 2 , k' 2 , and k 4 . 
Tables HT1 and Hill summarize our estimates. 
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To measure the k's, for a given value of q, for each 
of several system sizes (side length L a power of two 



FIG. 2: Winding angle variance for the largest contour when 
q — 1 as a function of the side length L of the box. Error bars 
on our estimates of the winding angle variance are shown, 
but are quite short and appear as points. A curve of the 
form K2 log L + a was least-squares fitted to this data and 
plotted here. The 95% confidence intervals (±1.96 standard 
deviations) for the parameters are K2 = 1.5002 ± 0.0023 and 
a = 0.55494 ± 0.013, consistent with k 2 = 3/2. The fit has a 
X 2 statistic of 23.65 with 23 degrees of freedom for a p-value 
of 0.42, so the fit passes the \ 2 test. In this case the fit is good 
all the way down to I/ m i n = 32. In some cases the fitted curve 
lies outside the 95% confidence interval at L m i n , indicating 
corrections to scaling, and in these cases we increased L m i n . 
These data are summarized in the first line of Table HT1 



We also conducted measurements for the minimum 
spanning tree (MST) with random edge weights. The 
paths of the MST are smoother and less windy than those 
of the UST (see Table EH . For MST it is unlikely that 
D f = l+Ki/8, so the MST path is not described by SLE. 

In conclusion, Eq. Q, which generalizes Duplantier's 
winding angle formula, is supported by both experiments 
and heuristic arguments. It would be interesting to see 
if Eq. {5} holds for SLE K . 

Acknowledgements: We thank Oded Schramm and 
Jane Kondev for valuable discussions. 
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TABLE II: Winding angle variance coefficient for the longest 
loop (K2), longest loop of the external perimeter {k! 2 ), and 
largest pinch of the l ong est loop (K4). When q = 4, we expect 
that log-corrections [43l |44T] affect the measured k' 2 and k' 4 . 
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